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We show using a realistic Hamiltonian-type model that definite outcomes of quantum measure- 
ments may emerge from quantum evolution of pure states, i.e quantum dynamics provides a de- 
terministic collapse of the wavefunction in a quantum measurement process. The relaxation of the 
wavefunction into a pointer state with classical properties is driven by the interaction with an envi- 
ronment. The destruction of superpositions, i.e. choosing a preferred attraction basin and thereby a 
preferred pointer state, is caused by a tiny nonlinearity in the macroscopic measurement apparatus. 
In more details, we numerically studied the many-body quantum dynamics of a closed Universe 
consisting of a system spin measured by a ferromagnet embedded in a spin-glass environment. The 
nonlinear term is the self-induced magnetic field of the ferromagnet. The statistics for the out- 
comes of this quantum measurement process depends on the size of the attraction basins in the 
measurement apparatus and are in accordance to Born's rule. 
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I. INTRODUCTION 

Quantum mechanics [1-3] is perhaps the most success- 
ful theory ever developed and provides very accurate de- 
scriptions of Nature. For example the gyromagnetic ra- 
tio is correctly predicted by quantum electrodynamics 
to twelve significant figures of accuracy [4]. There are 
at present no experimental findings that contradict the 
predictions made by the quantum theory! In spite of 
this enormous predicting power, the orthodox quantum 
theory [1-3] can only provide the probability for the out- 
comes of measurements. It can not describe how the 
complex quantum wavefunction describing the state of 
the system actually evolves into a single classical out- 
come, found in all measurement processes. To overcome 
this problem, the orthodox quantum theory postulates 
the discontinuous collapse of the wavefunction upon mea- 
surement, linking unitary quantum evolution to the clas- 
sical percieved world of single well defined events. This 
probabilistic limitation of the quantum theory has, since 
its conception, troubled many physicists [5-13]. For ex- 
ample Einstein stated in a letter to Born in 1926: "I, at 
any rate, am convinced that He is not playing at dice. " 
[14]. 

Essential insight to the measurement problem has been 
gained through the decoherence program [15-17]. Most 
importantly, realistic macroscopic systems are never iso- 
lated, but interacting with their environments. The inter- 
action with the environments forces the quantum wave- 
function of macroscopic systems to evolve into classical- 
like pointer states [15-17]. Decoherence theory explains 
the emergence of classical properties of systems coupled 
to their surrounding environments. But the theory has 
nothing to say about single outcomes! Decoherence al- 
ways refer to an ensemble average where the environmen- 
tal degrees of freedom are traced out, resulting in a re- 
duced density matrix stating the probabilities to find the 
system in different pointer states. There are no mech- 
anism in the Decoherence theory that chooses a single 



outcome among all the possible pointer states. To many 
practitioners of quantum mechanics this is not a prob- 
lem, since they interprete the wavefunction to only de- 
scribe ensemble averages. However, we here choose to in- 
terprete the wavefunction as representing an individual 
system. With this interpretation a state in superposi- 
tion of outcomes a and b cannot suddenly evolve to an 
outcome of a or b, even if decoherence is involved. 

In this work, we show from the quantum evolution of 
pure states that single outcomes of measurement events 
are possible as a consequence of a phase transition in the 
measurement apparatus. The wavefunction of the mea- 
surement apparatus is forced into classical-like pointer 
states by the interaction with the environment. The 
destruction of superpositions and, thereby, selection of 
a preferred pointer state is caused by an infinitesimal 
non-linearity in the macroscopic measurement apparatus. 
Which pointer state the wavefunction actually collapses 
into depends on fine details of the initial conditions. The 
probability for the wavefunction to collapse into a partic- 
ular pointer state is shown to be in close agreement with 
Born's rule. The statistics for the outcomes depends on 
the size of the attraction basins in the measurement ap- 
paratus which again is affected by the initial state of the 
measured object. 

To be more specific about our physical picture, con- 
sider first the quantum dynamics of a closed Universe 
consisting of a linear ferromagnet (A) embedded in a 
large environment (E) . Assume that the ferromagnet has 
an easy z-axis. Let the initial state of the linear ferromag- 
net be a superposition between a state where nearly all 
spins are pointing up and a state where nearly all spins 
are pointing down, \Ao) = \A^) + \Ay), and let the 
initial state of the environment be \Eq). When the ferro- 
magnet interacts with the environment, the ferromagnet 
will transfer its entropy and energy into the environment 
and relax towards its two ground states, \A^) and {A^), 
which turn out to be pointer states. But since the time 
evolution is linear the dynamics of the distinct branches 
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in the superposition is completely independent of each 
other. The wavefunction of the whole system then be- 
comes l^*) = + after some transient 
relaxation time. Already here we see that interactions 
with an environment do not provide single outcomes. 
Quantum superpositions survive the interaction with en- 
vironments. The closed Universe evolves into indepen- 
dent branches [9], here, a superposition of two pointer 
states. There are no mechanisms that select a preferred 
evolution branch and the ferromagnet remains in a su- 
perposition with zero average magnetization. We clearly 
sec the incapability of linear quantum mechanics to de- 
scribe everyday observations. Macroscopic ferromagnets 
are not in superpositions, their magnetization have pre- 
ferred directions. We show below that a small non-linear 
term in the ferromagnet may destroy the superpositions 
and select a preferred direction. 

Nonlinearity has been argued to exist intrinsically 
in quantum mechanics [18], for example at the Planck 
scale [19]. More importantly, macroscopic quantum co- 
herent systems are often effectively described by nonlin- 
ear Schrodinger equations [20]. For example the dynam- 
ics of the quantum coherent superconducting state is gov- 
erned by the nonlinear Ginzburg-Landau equation [21]. 
Another example is Bose-Einstein condensates where the 
mesoscopic quantum coherent wavefunction is governed 
by the nonlinear Gross-Pitaevskii equation [22]. The 
fractional quantum Hall state may also be described by 
an effective nonlinear equation [23]. From the "emer- 
gence" point of view, these kind of nonlinearity may nat- 
urally emerge as a consequence of the enormous number 
of interacting particles involved in the macroscopic co- 
herent state which may have less symmetry compared 
to the underlying many-body Hamiltonian [24]. Finally, 
nonlinearity may also rise naturally as a consequence of 
the interaction between two fluctuating quantum fields in 
the mean-field limit, see Appendix A. In this paper we 
view nonlinearity as modeling physical reality, not only 
a convenient (and necessary) approximation. 

A tiny nonlinearity will, to be shown below, force a 
macroscopic ferromagnet embedded in an environment 
to chose a preferred direction, consistent with experi- 
mental observations. For ferromagnets, one may argue 
for a physically sound, nonlinear term by considering the 
self-induced magnetic field created by the magnetization. 
Assuming that the spins in a ferromagnet are carried by 
charged particles, then in connection to the spin each par- 
ticle also possess a small magnetic moment, g^sS, where 
S is a dimensionless spin operator, g is the gyro magnetic 
factor and /xb is the Bohr magneton. The magnetic mo- 
ments create a magnetic field which may be, crudely, ap- 
proximated to B = /io5/^B(*| J2ieA where is the 
magnetic constant and j^) is the conventional normalized 
wavchmction of the system. Here, we have completely 
neglected spatial variations, fluctuations and other inter- 
nal degrees of freedom in the magnetic field. We believe 
that those additional degrees of freedom in B does not 
fundamentally modify the physical picture presented be- 



low. It is important to note that we use the operator 
(^'|...|\I') to denote a (inner) scalar product, no ensem- 
ble average is included in the process. The self-induced 
magnetic field may again interact with the spins in the 
ferromagnet giving rise to a nonlinear term in the Hamil- 
tonian which may, at the mean field level, be expressed 
as 

ifs = -M(*|^Sil*)-5^S„ (1) 

where the parameter jj controls the nonlinearity and has 
the dimension of energy. This non-linear term will favor 
the pointer state where the spins are parallel to the mag- 
netic field and disfavor all other pointer states. Here is 
the physical picture. Each measurement apparatus has a 
set of attraction basins, related to the set of eigenvalues. 
Each attraction basin has a stable fixed point that is the 
pointer state. Interactions with environments will force 
all the parts of a superposition that reside in an attrac- 
tion basin to relax into the basin's stable fixed point, i.e. 
the basin's pointer state. A self-induced magnetic field, 
B, above a certain threshold will single out a preferred at- 
traction basin and transform all other initially attraction 
basins into repulsive basins and their stable fixed points 
into unstable fixed points. Consequently, the measure- 
ment apparatus has now a single unique attraction basin 
with its corresponding pointer state. Hence, the ferro- 
magnet will, in the presence of a self-induced magnetic 
field together with an environment, evolve into a single 
unique pointer state. 

In the case of no coupling to an environment, there is 
no mechanism for the spins in the ferromagnet to relax 
and the ferromagnet is then forced to remain in a super- 
position regardless of nonlinearity. With an environment, 
however, non-zero fluctuations in the magnetization will 
create a magnetic field which, when strong enough, en- 
hances parallel spins and thereby increase the magneti- 
zation along that direction. Thus, environment induced 
fiuctuations may starts a self-enhancing process which fa- 
vors the pointer state with spins parallel to the direction 
of the initial fiuctuation. This process forces the sys- 
tem to choose a preferred attraction basin. During this 
process, the interactions with the environment also cause 
the wavefunction to fall into the fixed point of the chosen 
attraction basin, i.e. its pointer state. The probability 
for the wavefunction to end up in a given pointer state 
is governed by the size of the attraction basins, which 
for a ferromagnet alone is of course 50% for ending as 
I*) = \A^)\E^) and 50% for ending as |*) = \Ai)\Ei). 
We will discuss the the size of the attraction basins in 
more details below. In summary, a non-linear term in the 
Hamiltonian as Eq. (1) may somewhere along the time 
evolution of the ferromagnet chose a preferred attraction 
basin. In the same time, interactions with the environ- 
ment will force the wavefunction to fall into an unique 
classical-like pointer state which is the stable fixed point 
of the chosen attraction basin. Which attraction basin 
that is actually chosen depend on small details of the ini- 
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tial state of the ferromagnet and the environment. An 
infinitesimal smaU change in the initial state may force 
the wavefunction to collapse in a entirely different pointer 
state. However, the dynamics is completely determinis- 
tic. Given the same initial condition, the wavefunction 
will always choose the same attraction basin and collapse 
into the same final pointer state. In this sense, quantum 
mechanics is a deterministic theory. 

Nonlinearity in the ferromagnet can in principle be in- 
finitesimally small and still do its job in the selection of 
the preferred attraction basin. The reason is as follows. 
For simplicity, consider a ferromagnet in 2 dimensions 
with the Hamiltonian 

Hj = -J^S^S^ (2) 

{i,3) 

where {i,j) denotes nearest neighboring spin-i spins. 
This ferromagnet has two degenerate ground states, | tt 
... t) and I U ... i) with the energies Ei; = = -JNa/2 
where Na is the number of spins in the ferromagnet. 
The ground states are separated by an energy barrier 
with height —JNa/2. A magnetic field given by Eq. (1) 
will lift the degeneracy between E^ and E^ by AE — 
El — El = ijlNa/^. Note that this argument also ap- 
plies for the case where the magnetic field decays a as a 
function of distance from its sources provided that the ex- 
change interaction between the spins decays much faster. 
Thus for large Na, the energy from the initially fluctu- 
ating magnetic field will not only define a global mini- 
mum, but also eventually grow strong enough in order to 
completely remove the local minimum connected to the 
anti-parallel ground state forcing both branches of the 
wavefunction to fall into the same fixed point defined by 
the global minimum. Therefore, an infinitesimally small 
nonlinear parameter /z is sufficient to define a global mini- 
mum and remove all other local minima for a macroscopic 
ferromagnet. 

We will from now on speak of the ferromagnet {A) as 
an measurement apparatus, referring to the system as our 
measurement object. Including a measurement object, 
that in our case is a system spin Sgys interacting with the 
ferromagnet, will alter the size of the attraction basins 
of the ferromagnet. Here we provide an estimate for the 
probabilities showing that this change of size provides 
statistics in accordance with Born's rule. First, let us ex- 
clusively focus on the case of a ferromagnet embedded in 
an environment with no measurement object. For clar- 
ity we again restrict the discussion to an Ising ferromag- 
net with spin-i spins. Define the dimensionless magnetic 

field for the ferromagnet as B{t) = {^{t)\ J2^<.A SH'^it)), 
and let the initial magnetization of the ferromagnet be 
zero, B{t=0) = 0. Interactions with environments force 
B{t) to deterministically fluctuate back and forth around 
zero. Assume that when B reaches a threshold, zkBc, the 
nonlinear term dominates over the fluctuations and starts 
a self enhancing process which ends up in a macroscopic 
stable magnetization, sec Fig. 1. A measurement appara- 
tus that can measure the direction of a single spin must 



be sensitive to a flip of a single spin. Thus, Be — ±i. In 
other words, when the environment has flipped one spin 
in the ferromagnet, nonlinearity will set in and force all 
other spins to align to that spin. Note that this value 
of Be applies only for apparatus that are able to detect 
the direction of a single spin-^ spin. If the initial states 
of the ferromagnet and the environment are symmetric 
with respect to the up and down pointer states, B{t) will 
fluctuate symmetrically around zero. Hence, the proba- 
bility for that B(t) reaches B^ first and the ferromagnet 
will eventually choose the up state is 0.5. Similarly, the 
probability for that B{t) reaches —Be first and the ferro- 
magnet will eventually end up in the down state is 0.5, 
see Fig. 1. Note that the dynamic from a state with zero 
magnetization into a state with macroscopic magnetiza- 
tion closely resembles the dynamics of a phase transi- 
tion. The state of the system evolves from a state with 
the same symmetry as the underlying Hamiltonian into 
a state with less (broken) symmetry [24] . 

Now, add a system spin to be measured by a ferromag- 
net embedded in an environment. Let the system spin be 
in the initial state 

\Ssys{t^O)) ^ a\ ^) + f3\ i) (3) 

where a and /? are complex scalars. Let for simplicity the 
initial state of the ferromagnet be in an antiferromagnetic 
configuration 

l^o) = I TiTi - Ti). (4) 

Of course, other states including those with superposi- 
tions which has zero magnetization will do the job. Af- 
ter a short time the system spin has interacted with the 
ferromagnet and flipped its spin leading to a combined 
state where one spin in the ferromagnet is fiipped, e.g 

* = a| i)| TTTi...Ti)+/?l T)l iiti-Ti). (5) 

The dimensionless magnetic field is now B = ^{a^ — 0^) 
so the starting point of B{t) is no longer at zero but at 
\{oi^ — (3"^), see Fig. 1. We see also that the size of the 
"up" attraction basin is | -|- ^(a^ — /3^) while the size 
of the "down" attraction basin is ^ — ^(a^ ^ Z?^)- The 
probability ratio for B to reach B^ first versus —Be is 
therefore, see Fig. 1, 

In the last equality we have used the normalization con- 
dition + = 1. Thus, the probability for the ferro- 
magnet to end with a stable up state is Pf = a^. And, 
the probability for the ferromagnet to end with a sta- 
ble down state is Pi = f?, which are in accordance with 
Born's rule. 

Our approach should not be confused with that of 

the dynamical reduction program [12, 25], where the 
Schrodinger equation is altered by inclusion of stochastic 
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time 



Bc{l)=-i- 



Figure 1: Schematic figure of the dimensionless magnetic field 

as a function of time. Due to the initial state of the system 
spin |5sys(t = 0)) = a| t) + P\ I), the dimensionless magnetic 
field start of its deterministic random walk from ^{a^ — P'^). 

When B{t) reaches one of the two threshold values, non-linear 
effect will dominate over the fluctuations and force the ferro- 
magnet to chose a direction. 



nonlinear terms in order to achieve dynamical or sponta- 
neous localisation of wave packets. The basic idea behind 
the dynamical reduction program is that the wavefunc- 
tion randomly collapses all the time, where the rate of 
collapse is related to the mass of the system. This applies 
to all particles whether isolated or interacting [12, 25]. In 
our approach, an isolated system remains unperturbed 
until the different systems start to interact. It is the 
asymmetry between the size of the environment and the 
rest that causes a classical appearance of the state of the 
measurement apparatus and the measured object upon 
measurements. 

Ref. [26] considered a self-collapse mechanism where 
the quantum wavefunction automatically exhibit dynam- 
ical collapse, without any measurement apparatus nor 
environment, due to nonlinearity originating from non- 
Abelian gauge fields. Our description of a measurement 
process is completely different from Ref. [26]. We find 
that the dynamics of the wavefunction strongly depend 
on the existence of a macroscopic measurement appara- 
tus embedded in an environment that is able to absorb 
entropy and energy from the measurement apparatus. 

Ref [27] studied a collapse of the quantum wavefunc- 
tion due to a spontaneous symmetry breaking process in 
the measurement apparatus, an anti-ferromagnet in the 
thermodynamical limit. To induce the collapse Ref [27] 
used a time dependent, fluctuating, symmetry breaking, 
non-unitary, staggered magnetic field with an undeter- 
mined origin. A major difference between our model and 
the model in Ref [27] is the importance of the environ- 
ment in the collapse of the wavefunction. We studied 
the evolution of a pure wavefunction describing a closed 
Universe. 

Note that the results of this work by no means contra- 
dict the predictions of quantum mechanics. Merely, the 



results shows that quantum mechanics may be regarded 
as a complete theory which can describe measurement 
processes and predicts single outcomes of classical alike 
states without the collapse postulate. In other words, the 
results show that quantum mechanics is a deterministic 
theory. The price to pay is a tiny non-linearity effectively 
only in systems of macroscopic sizes. 



II. MODEL 

To support the ideas discussed above in the introduc- 
tion, we perform numerical modelling of an idealized, but 
realistic, model of a quantum measurement process. Our 
model is composed of a system spin, Ssys, being mea- 
sured by a ferromagnet (A) embedded in a spin-glass en- 
vironment {E). The system spin is a spin-i spin. The 
measurement apparatus is a ferromagnet consisting of a 
number = 4 or = 8 spin-i spins. In addition to 
the traditional spin-spin exchange couplings, the spins in 
the ferromagnet also interact with its self-induced mag- 
netic field, giving rise to a weak non-linear term. Finally, 
the system spin and the ferromagnet are embedded in 
an environment consist of a Ne = 15 spin-i spins with 
random, frustrated, spin-glass interactions. 



A. Hamiltonian 

Our model Hamiltonian may be written as 

H = Ha + He + Hae + Hsa+Hse + Hb (7) 

rr \ ^ \ ^ ra oq qcx. 

Ha - - 2_^Jijbi Sj 

ieAjeE a 
iJeE a 

\ " \ " oa oa ja 

/ J / J ^sys i 



Hae 
He 

HsE 

HsA = ^TiS^ygS, 
Hb 



ieE a 



ieA 



= -i,{^\Y,sn^) E^i 



ieA 



jeA 



where where S"", /" and Sfy^ are dimensionless spin-i op- 
erators in the ferromagnetic measurement apparatus, the 
spin-glass environment and the system spin to be mea- 
sured, respectively. In Eq. (7), the index a = x,y,z runs 
over the three components of the spin operators. The ex- 
change couplings Jj" , Qfj and control the interaction 
between spins in the ferromagnet, the spin-glass envi- 
ronment and between the spins in the ferromagnet and 
the environment, respectively. Furthermore, Fj controls 
the interaction between the system spin and the appara- 
tus while Of controls the interaction between the system 
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spin and the environment. In more physical terms, the 
measurement apparatus is modelled by a ferromagnet, 
Ha- a spin glass, He, serves the role of an environment 
which absorbs the energy and entropy of the ferromag- 
net allowing it to relax towards its ferromagnetic ground 
states. The systcmi to be measured is a single spin-i spin, 
Ssys, which role is to tilt the up-down symmetry of the 
ferromagnet. The interaction between the ferromagnet 
and the spin glass is modelled by Hae- While Hsa de- 
scribes the interaction between the system spin and the 
ferromagnet. To make the measurement process more re- 
alistic, wc have also a coupling between the system spin 
and the environment, Hse- 

The Hamiltonian given by Eq. (7) without the non- 
linear term, Hb, is linear and will not be able to conduct 
any measurement process as discussed in the introduc- 
tion and shown in more details below. This linear part of 
the Hamiltonian can not choose a preferred pointer states 
among the possible pointer states for the ferromagnet. To 
model a ferromagnet that shows the behavior found in ex- 
periments, we introduce a non-linear term, Hb , which de- 
scribes the interaction between the self-induced magnetic 
field and the spins in the ferromagnet. Physically, the 
self-induced magnetic field (B) originates from the mag- 
netization of the ferromagnet, B oc M oc X]ieA(*l'5'l'l*)- 
Again, |^) is the normalized wavcfunction for the whole 
system, i.e. closed Universe. We use the parameter 
to control the non-linear self coupling between the mag- 
netic field and the spins Si in the ferromagnet. In the 
macroscopic limit, an infinitesimal small ji will force the 
ferromagnet to choose a direction. Note that these kind 
of one-particle non-linear interaction can also be argued 
to originate from linear terms like Si-Sj treated in a mean 
field level Si-Sj {Si) ■ Sj. This kind of mean field 
approximation has been used with tremendous success 
in other macroscopic many-body quantum systems as 
superconductivity [21], superfluidity [22] and fractional 
quantum hall states [28-30]. There, the mean field so- 
lution accurately describes the properties of complicated 
many-body quantum wavefunctions in the macroscopic 
limit. Another source for the nonlinear term is the in- 
teraction between and a magnetic field treated in the 
mean-field approximation, see Appendix A. Note that 
the Schrodinger equation used, Eq. 7, is of Hamiltonian 
type even in the presence of a nonlinear term. The energy 
is conserved during time evolution, see Appendix A. 

Now, to the parameters of the Hamiltonian, Eq. 7. We 
assume J?- = J for nearest neighbors and J/\/2 for next 
nearest neighbors. The spin configuration of the ferro- 
magnetic measurement apparatus is rectangular (cubic) 
when composed of 4 (8) spins, respectively. For clarity wc 
will, henceforth, use J as our unit of energy, i.e. J = 1. 
All other energies are measured with respect to J. We use 
anisotropic coupling Jf^ = Jfj = ^J^^ =7; where 7 < 1, 
to make the ferromagnet Ising-like. This is in order to 
reduce the infinite set of pointer states for a Heisenberg 
ferromagnet to a set of two pointer states, up and down, 
for a more transparent physical picture. The exchange 



couplings Vtf, , A°^- and 0f are uniform random numbers 
in the range fi], [-A, A] and [-9, 6], for ah a. The 
coupling between the system spin to be measured and the 
measurement apparatus = 1, for all i E A. 

A similar Hamiltonian without the system spin and the 
non-linear self induced magnetic field B has been stud- 
ied in Rcf. [31]. There, they study, numerically, how the 
spin glass (E) relaxes the ferromagnet (A) [31-33]. It was 
found that frustrated spin-glasses where very effective 
with respect to relaxation and decoherence of the cen- 
tral system, even for small numbers of spins {Ne ~ 10) 
in the environments. Spin glasses are therefore ideal as 
environments for small spin systems. 

B. Preparation and numerical method 

This section describes the preparation of the initial 
state and the procedure for the numerical time integra- 
tion of the Schrodinger equation. 

First, we prepare an initial state 

= 0)) = \Ssys{t = 0)) ® \Ao) ® \Eo). (8) 

The system spin to be measured is prepared in a general 
superposition 

\S,ys{t = 0)) = \ T) cos^ + 1 i) sin^e'^, (9) 

where the phase angle (p G [0, 27r] is not so important for 
the physical picture and henceforth set to 0. The ferro- 
magnet is for simplicity prepared in an antiferromagnetic 
state \Ao) = | fj. ... fj,). Of course, other states with zero 
magnetization would do the same job. The environment 
is prepared in a state close to the spin-glass ground state 
\Eo) computed by the Lanczos method [34]. 

Next, we iteratively time integrate the Schrodinger 
equation = H'^ using the finite difference method, 

^{t + dt) =^{t)-iH^dt, (10) 

starting from the initial state, Eq. (8). The Hamiltonian 
H is given by Eq. (7). For clarity, we have set h = 1 in 
Eq (10) and throughout the article. Each iterative step 
in Eq. (10) is found numerically by the Chebyshev expan- 
sion method [35-37]. The method is based on a polyno- 
mial expansion of the propagator U{t,to) = e'(*~*°)^ w 

K 

Jo{t) + Jk{t)Tk{H), where Tk{x) is the Chebyshev 

k=l 

polynomial of the first kind and Jk {x) is the Bessel func- 
tion of integer order k. The Chebyshev moments Tk{H) 
is computed by the recursive application of the Hamil- 
tonian operator H on the wavcfunction l^*). We have 
verified that the numerical method preserves the norm 
of the wavefunction, {'9{t)\^{t)) for all t. Furthermore, 
the conserved energy of the closed universe 

Eu = {^\H - ^Hb\^} (11) 
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has also been verified to be conserved in the numerical 
calculations. A short derivation of this form of the con- 
served energy for our model Eij is given in Appendix A. 

The Hilbert space of the composite system H = Hs (S) 
Ha® He is of dimension 2^, where L is the total number 
of spins, system spin plus ferromagnet plus spin-glass. 
Present computing power limits us to model systems 
where L < 24. 

The following parameters are fixed in all simulations: 
A = 0.3, Vt = 0.8, e = 0.5, Ne = 15. The remaining 
parameters: Na — 4: ^ ^ = 12.0 or Na = 8 ^ ^ — 6.0. 
The nonlinear parameter /i is varied such that the effect 
of maximum magnetization on each spin in the ferromag- 
net is constant, (iNa = 48.0. This allows direct compar- 
isons between the measurement results from ferromag- 
nets of different sizes. Since our closed Universe is very 
small, the non-linear parameter had been set relatively 
large to optimize the measurement properties of our tiny 
measurement apparatus. In the macroscopic limit the 
non-linear parameter may be infinitesimal small and still 
do the same job, as argued for in the introduction. 



III. RESULTS 

In this section we present the numerical results for 
the model given by Eq. (7). The results support our 
claim that the wavefunction of the measurement appara- 
tus, in the presence of a small nonlinear term, is forced 
to choose one of the two classical pointer states \A^) or 
l^j^). Without the non-linear term, the ferromagnet re- 
mains in the superposition regardless of the interaction 
with the environment. We also present statistics for mea- 
surements with different initial configurations for various 
6, and show that the presence of a measurement object 
will alter the size of the attraction basins in the mea- 
surement apparatus in a way such that the outcome is 
consistent with Born's rule to the numerical precision. 



A. Time evolution 

It is difficult to extract useful information directly from 
the complex many-body quantum state of the whole sys- 
tem which is, here, a 2^° or 2^^ dimensional complex 
vector. In order to characterize the ferromagnet, we use 
the dimensionless magnetization 

Mit) = J2{^{t)\Sn^{t)) (12) 

which determines the degree of ferromagnetic order in 
the measurement apparatus. In addition, we use the ex- 
change energy 

E{t) = - J2 J^,{'^it)\S^■S,\^it)) (13) 

i,jl£A 



to characterize the relaxation of the ferromagnet. To 
characterize the state of the system spin we use the ex- 
pectation value 

(Slysit)) - (*(i)|5f^Jvi,(t)). (14) 

Again, note that we use the operator to denote a 

(inner) scalar product, no ensemble average is performed. 
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Figure 2; Typical time evolution of the magnetization Al 
(solid) and the exchange energy E (dashed) for a linear ferro- 
magnet and the expectation value of the system spin (dashed 
dot). The initial state of the system spin is {Saysit — 0)) = 
'^il T) + I i))- The number of spins in the ferromagnet is 

Fig. 2 shows a typical time evolution of the magne- 
tization AI and the exchange energy E characterizing 
the state of a linear ferromagnet, i.e. ^ = 0. The sys- 
tem spin is initially prepared in the state \Ssys{t = 0)) = 
"^(It) + li))- We see that the exchange energy of the 
ferromagnet decreases rapidly, as the ferromagnet relaxes 
from the anti-ferromagnetic state towards its ferromag- 
netic ground states, \Af) and |^j,). Note that due to the 
limited number of spins used in the ferromagnet, Na = 4, 
the interaction with the large environment prevents the 
ferromagnet to relax all the way down to the ferromag- 
netic ground states where E = —1.35. We also see in 
Fig. 2 that the magnetization M shows small fluctua- 
tions, but no sustained net magnetization characterizing 
the expected classical ground state. Fig. 2 also shows 
that the expectation value of the measured system spin 
is zero, indicating that jS'sj^s) remains in a superposi- 
tion between ||) and |J,). A more careful examination 
of the stationary state shows that the composite wave- 
function of our closed universe is still in a superposition. 
Thus, the magnetization fluctuates around zero and the 
ferromagnet does not choose any preferred pointer state. 
From Fig. 2, one may conclude that linear ferromagnets 
are not able to conduct any measurement. The super- 
position persists regardless of the interaction with the 
environment. The time evolution 

(I ^)+\i))\Ao)\Eo) ^ I m^m)+\mim), (15) 

describes the dynamics in the linear case. The initial 
state of the system becomes entangled with the appara- 
tus and the environment. The two branches of Eq. (15) 
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evolve completely independent of each other, and the ini- 
tial superposition has through interaction with its envi- 
ronment extended to our entire closed universe. Thus, 
decoherence, in the sense of relaxation to the pointer 
states, alone is not sufficient to remove the initial su- 
perposition. It should be noted that the persistence of 
the superposition in this case is entirely due to the exact 
linearity of the dynamical equation, and has nothing to 
do with the size of the apparatus. It is straight forward 
to show that superpositions survive any linear finite sized 
ferromagnet, see Appendix B. 
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Figure 3: A typical time evolution of the magnetization M 
(solid) and the exchange energy E (dashed) of a non-Unear 
ferromagnet, and the expectation value (Ssys) (dashed dot). 
The system spin is initially prepared in the state \Ssy3{t — 
0)) = -1=(| T) + \ i)). The number of spins in the ferromagnet 
is Na = 4. 
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Figure 4: Same as in Fig. 3. The only difference is the initial 
state of the environment. 

Typical time evolutions of the magnetization M of a 
non-linear ferromagnet is shown in Figs. 3 and 4. After 
initial fiuctuations due to the environmental interaction, 
one of the two possible magnetization directions is en- 
hanced due to the feedback provided by interaction with 
the self-induced magnetic field B. A "phase transition" 
in the ferromagnet takes place, where the ferromagnet 
relaxes into the pointer state of the attraction basin se- 



lected by the emerged magnetic field. We have quoted 
the term "phase transition" here since our tiny ferromag- 
ncts consist of only 4 spins and phase transitions happens 
in principle only in the thermodynamical limit. Rather, 
finite size effects [38] are expected to dominate the behav- 
ior of our small ferromagnets. This can be seen in the 
actual values for M{t), E{t) and Ssys{t)- They are all 
significantly less than the values expected in the thermo- 
dynamic limit which are: M = Q.^-^Na — 2.0, E — —1.35 
and Ssys = 0.5. However, the stationary values of M, E 
and Ssys are stabilized well beyond the noise and one may 
clearly conclude that a measurement has been conducted 
by an admittedly rather poor measurement apparatus. 

In more physical terms, when the magnetization 
reaches a critical value, the non-linear term sets in and 
turns all attraction basins, except the chosen one, into 
repulsive basins. From this point there is no turning 
back, all superpositions in the wavcfunction of the closed 
Universe are forced to evolve such that the state of the 
ferromagnet becomes the unique pointer state chosen by 
the magnetic field. Thus, as a result of the measure- 
ment dynamics the spin expectation value of our system 
{Ssys) will follow the direction of the established mag- 
netization M . The state of the measurement object has 
been driven into the state corresponding to the final out- 
come of the measurement apparatus. The final result is 
in correspondence to the result obtained by the collapse 
postulate of orthodox quantum theory where an initial 
state (I t) + l i))® l^o) is instantaneously and discontinu- 
ously collapsed upon measurement to either | ]) ® or 
I I) ® \Ai). We might therefore conclude that a quantum 
measurement of a single event has been accomplished. 
The total wavefunction is 

\^)^\])\A^)\E^) (16) 

after a measurement where the up pointer state is cho- 
sen. Note that the final state of the environment 
still contains superpositions inherited from the initial su- 
perposition of the measured spin. 

The only difference in the initial state \^) of the com- 
posite system between Fig. 3 and Fig. 4 is the different 
initial state of the environment. Note that the different 
initial states of the environment are completely equiva- 
lent, they all belong to the set of degenerate ground states 
of the spin glass. We see that a small difference in the 
initial state of the environment does in this case alter the 
dynamics of the measurement process completely. In this 
example the slight altering of the initial state resulted in 
a different final pointer state of the measurement appa- 
ratus. It is this sensitivity to the initial state that gives 
rise to the indoctrinated randomness in orthodox quan- 
tum mechanics, and to Born's rule. 

B. Statistics 

According to Born's rule [39], in an ideal quantum me- 
chanical measurement the probability for our measure- 
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ment apparatus to end in the state of positive magne- 
tization |v4|), when the object is prepared in the state 
\S,y,it^O)) of Eq. (9) is 

P|(6i) = cos^ei. (17) 

Correspondingly, the probabihty for ending in the nega- 
tive magnetization state is 

Pl{e)=sin^9. (18) 

To obtain measurement statistics, we run a number of 96 
independent simulations for each choice of 9. For each 
simulation, a different spin glass ground state is used as 
the initial state for the environment. 
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Figure 5: Probability for a ferromagnet of 4 spins to end in 
an up pointer state as a function of the angle 9 of the initial 
system spin state |S'sj,s(t = 0)) = | |)cosS+ | J,)sin0. Solid 
line shows ideal statistics according to Bom's rule. 

The probability for the measurement apparatus and 
the system spin to end up as a function of 6 for a fer- 
romagnet consist of Na = 4 spins is shown in Fig. 5. 
Here, the non-linear parameter is fi = 12.0. We have 
verified that Pi{9) + Pii9) = 1. We see in Fig. 5 that 
the probabilities of obtaining the pointer state | t)I^T) 
the outcome of the measurement process resemble Born's 
rule. The discrepancies at small and large 9 are due to 
the huge finite size effect in the tiny measurement appa- 
ratus. The subsystem of system spin plus apparatus are 
so small, that the interaction with the large environment 
affects the measurement result. Needless to say, a 4-spins 
ferromagnet is far from a perfect measurement appara- 
tus. It is, for us, surprising that this tiny ferromagnet 
actually provides results so close to Born's rule. One 
reason may be due to the up-down symmetry of the fer- 
romagnet and of the environment. Since the subsystem 
of system spin plus ferromagnet is so small, the "noise" 
from the environment dominates and the probability for 
ending up is close to 50% regardless of the initial state of 
the system spin. Thus, P^{9 — 45) = 0.5 is guaranteed 
by symmetry. 

To show that larger ferromagnets will actually provide 
better measurement results, we have also carried out sim- 
ulations with a system with a Na = 8 spins ferromag- 
net. Here, the non-linear parameter is fi = 6.0. Recall 
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Figure 6: Probability for a ferromagnet of 8 spins to end in 
an up pointer state as a function of the angle 6 of the initial 
system spin state |5'3ys(i = 0)) = | t)cos0-|- | J.)sin6. Solid 
line shows ideal statistics according to Born's rule. 



that, this value of fj, is chosen to conserve the effect of 
maximum magnetization on each spin in the ferromag- 
net, (J-Na = 48, such that a direct comparison with the 
Na = 4 case is possible. All other parameters are iden- 
tical as above. As seen in Fig. 6, the probability for the 
ferromagnet to end in an up state follows Born's rule 
more closely compared to the Na — 4 spins ferromagnet. 
This is expected, since the finite size effect is smaller 
here. By comparing the cases Na — A ^ fi — 12.0 and 
Na = 8 ^ — 6.0, we see that the nonlinear param- 
eter fi decreases for increasing size of the ferromagnet. 
Thus, for a macroscopic system the non-linear parame- 
ter fj, may be infinitesimally small and still do its job: 
selecting a preferred attraction basin and consequently a 
pointer state. 

Obviously the obtained statistics is highly dependent 
on the choice of parameters, for example on the strength 
of the nonlinear interaction /i. We have chosen the pa- 
rameters of our model in such a way that our mea- 
surement apparatus is optimized for measuring the z- 
component of a single spin. 



IV. CONCLUSION 

We show from a model study of a simple closed Uni- 
verse that definite outcomes of quantum measurements 
can emerge continuously from pure quantum evolution. 
Each measurement apparatus has a set of attraction 
basins corresponding to the set of eigenvalues. The se- 
lection of a preferred attraction basin in a measurement 
process is caused by an infinitesimally small nonlinearity 
in the measurement apparatus. Interaction with the en- 
vironment then forces all superpositions of the measured 
object and the measurement apparatus to fall into the 
same fixed point (pointer state) of the preferred attrac- 
tion basin. The dynamics of this measurement process 
which strongly resembles the dynamics through phase 
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transitions is entirely deterministic! Given the same ini- 
tial condition the wavefunction will always fall into the 
same pointer state. However, small changes in the initial 
state, for example in the environment, may cause the 
wavefunction to fall into a completely different pointer 
state. Thus, the uncontrollable degrees of freedom in the 
environments account for the intrinsic statistical behav- 
ior of the quantum measurement process. The probabil- 
ity for falling into a certain pointer state is govern by the 
size of the attraction basins of the measurement appara- 
tus and is shown to be in close agreement with Born's 
rule. 

How the dynamics of non-local experiments, e.g. the 
EPR experiment [40] , will be affected by non- linear mea- 
surement apparatus is an interesting future study. 

To conclude, we have presented a model where the 
probabilistic behaviour of Quantum mechanics has ex- 
actly the same origin as the random outcomes of a rolling 
dice, i.e. its sensitive behaviour of uncontrollable initial 
conditions. Thus, even though neither He nor Quantum 
Mechanics are playing at dice, we physicists still have to. 
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Appendix A: CONSERVATION OF ENERGY 

In this appendix, we derive the expression for the con- 
served energy for our model, Eq. (11). The non-linear 
term is shown to introduce an extra term in the expres- 
sion for the conserved energy in addition to the standard 
average of the Hamiltonian. The appendix also shows 
how a non-linear term as the one used in our model may 
naturally arise from the interaction between two different 
fields. 

For clarity, we consider a simple quantum field model 
having the basic properties of our model, Eq. (7). The 
Lagrangian of our quantum field model with Einstein 
sum-convention over repeated indices assumed can be 
written as 

L = - *:/f„6*6 - ^ocKKb^b + - V{^). 

(Al) 

Here, is a complex field representing the wavefunction 
and Hab is the Hamiltonian for "J* where a and h are multi- 
indices which may include spin number, spin direction 
and also continuous variables like space coordinated etc. 
In Eq. (Al), <I> is a scalar field which may be associated 
with for example a harmonic oscillator with mass e or in 



certain limit the electromagnetic field. The multi-index 
a denoting all the variables necessary to describe the field 
$ which lives in a potential V($). 
Prom the Euler-Lagrange equation 

we deduce the dynamical equations 

= + $,r«(,*6 (A4) 
'^" + a¥;^^*) = (A5) 

Without the coupling between the fields when F";, = 0, 
Eq. (A4) becomes the conventional Schrodinger equa- 
tion for \1/ in the first quantized language. Furthermore, 
Eq. (A5) take the form of Newton's second law for the 
field 

To deduce the conserved energy we now assume that 
the field $ is, similar to the electromagnetic field, mass- 
less i.e e = and that the potential is of the form 

F($) = (A6) 

corresponding to the potential energy for the electromag- 
netic field ^-B^. Equation (A5) is then simplified to the 
constraint equation 

= -**r«6*6, (A7) 

which inserted into the Lagrangian Eq. (Al) give the 
effective Lagrangian for the field 

L' = Z^li>a - KHab^b - i$a*:r«(,*6. (A8) 

We can now deduce the conserved energy 

E = KHab'i^b + l^aK^'^b'^b- (A9) 

The conserved energy in our model, Eq. (7), is then 

Eu = - iHsl'i^) (AlO) 

if one relate the interaction term —^a^^T"f^^b 
in Eq. (Al) with our non-linear term Hb = 
-m(X; (<??)) E and $ with J^i^i)- Thus, we see 

that the nonlinear term Hq gives rise to an additional 

term —^Hg in the expression for the conserved energy 
compares to the conventional E = (H). 

More importantly, one see that the interaction between 
the fields $ and ^! give rise to a nonlinear term for the 
dynamic of in the mean field limit, i.e when $ ($). 
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Appendix B: LINEAR TIME EVOLUTION 

We will in this section show that superpositions survive 

interaction with arbitrary macroscopic measurement ap- 
paratus as long as the time evolution is linear. Consider 
a system to be measured, initially prepared in the super- 
position If) + II). The system is being measured by a 
measurement apparatus A, interacting with an environ- 
ment E. We can write the initial state of the composite 
system as 

\^) = {n) + \i))^\Ao)®\Eo) (Bl) 

According to linear quantum mechanics, time evolu- 
tion is governed by the time evolution operator U{t,to), 
acting linearly on any quantum state 

U{t, fo)(*A + *b) = U{t, io)*A + U{t, io)*B. (B2) 



Thus the time evolution of branch A of the wavefunction 
is completely independent of branch B. 

When applied to a quantum measurement linear quan- 
tum evolution lead to the von Neumann chain of infinite 
regress [3] 

(IT) + li)) j^o) \Eo) - IT) lA) \E^) + IT) l^i) \E^) . 

(B3) 

Since the environment (E) can in principle denote every 
degree of freedom in the universe, including observers, 
we see that linear time evolution result in the initial su- 
perposition of the system being measured eventually ex- 
tending to the entire universe. This time evolution (i.e 
linear quantum mechanics without the collapse postu- 
late), naturally lead to the Everett relative state inter- 
pretation [9]. Where the universe continuously split in 
distinct independent branches. 
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